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Abstract 

We consider spatially homogeneous marked point patterns in an unboundedly expand- 
ing convex sampling window. Our main objective is to identify the distribution of the 
typical mark by constructing an asymptotic x 2 -goodness-of-fit test. The corresponding 
test statistic is based on a natural empirical version of the Palm mark distribution and 
a smoothed covariance estimator which turns out to be mean-square consistent. Our 
approach does not require independent marks and allows dependences between the mark 
field and the point pattern. Instead we impose a suitable /3-mixing condition on the un- 
derlying stationary marked point process which can be checked for a number of Poisson- 
based models and, in particular, in the case of geostatistical marking. Our method 
needs a central limit theorem for /3-mixing random fields which is proved by extending 
Bernstein's blocking technique to non-cubic index sets and seems to be of interest in its 
own right. By large-scale model-based simulations the performance of our test is studied 
in dependence of the model parameters which determine the range of spatial correlations. 
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1 Introduction 

Marked point processes (MPPs) are versatile models for the statistical analysis of data 
recorded at irregularly scattered locations. The simplest marking scenario is independent 
marking, where marks are given by a sequence of independent and identically distributed 
random elements, which is also independent of the underlying point pattern of locations. A 
more complex class of models considers a so-called geostatistical marking, where the marks 
are determined by the values of a random field at the given locations. Although the random 
field usually exhibits intrinsic spatial correlations, it is assumed to be independent of the 
location point process (PP). However, in many real datasets interactions between locations 
and marks occur. Moreover, many marked point patterns arising in models from stochastic 
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geometry such as edge centers in (anisotropic) Voronoi-tessellations marked by orientation or 
PPs marked by nearest-neighbour distances do not fit the setting of geostatistical marking. 
Statistical tests for independence between marks and points are e.g. discussed in [8, 9, 23, 25]. 
A frequent approach to investigate dependences in marked point patterns is based on mark 
variogram and mark covariance functions. Recently, asymptotic normality of empirical ver- 
sions of these functions with applications to mark correlation analysis has been studied in 
[10, 11, 13]. The main goal of this paper is to investigate estimators of the Palm mark 
distribution P M in point patterns exhibiting correlations between different marks as well as 
between marks and locations. The probability measure P M can be interpreted as the distri- 
bution of the typical mark which denotes the mark of a randomly chosen point of the pattern. 
For any mark set C we consider the scaled deviations Z k {C) = y/\W k \{ (P M ) k ( C ) ~ P m( C ) ) 
as measure of the distance between P M and an empirical Palm mark distribution {P°[)k ■ 
Under appropriate strong mixing conditions we are able to prove asymptotic normality of the 
scaled deviation vector = (Z k (C\), . . . ,Zk(Ce)) T when the observation window W k with 
volume \Wk\ grows unboundedly in all directions as k — > oo. The proof relies on Bernstein's 
blocking method, see e.g. [4, 21], which so far has been applied only to sequences of cubic 
or cubelike windows W k , see e.g. [11, 12]. By means of some convex-geometric arguments it 
turns out that the blocking method is indeed applicable to any increasing sequence of convex 
observation windows W k with unboundedly growing inball radii. In addition we discuss con- 
sistent estimators for the covariance matrix of the Gaussian limit of Z&. This enables us to 
construct asymptotic x 2 -goodness-of-fit tests for the Palm mark distribution P M . By means 
of computer simulations we study the convergence of first and second type errors of the tests 
for growing observation windows in relation to the range of dependence of the MPP. In this 
way we demonstrate the practicability of the tests in analysis of real data. A promising field 
of application of our testing methodology could be the directional analysis of random sur- 
faces. Based on our results one can e.g. consider Cox processes on the boundary of Boolean 
models, marking them with the local outer normal direction and testing for a hypothetical 
directional distribution. This allows to identify the rose of directions of the surface process 
associated with the Boolean model and represents an alternative to a Monte-Carlo test for 
the rose of direction suggested in [2]. The occurring marked point patterns differ basically 
from the setting of independent and geostatistical marking, for which functional central limit 
theorems (CLTs) and corresponding tests have been derived in [16, 22]. Our paper is orga- 
nized as follows. Section 2 introduces basic notation and definitions. In Section 3 we present 
our main results, which are proved in Section 4. In Section 5 we briefly discuss some models 
satisfying the assumptions needed to prove our asymptotic results. In the final Section 6 we 
study the performance of the proposed tests by large-scale simulations. 



2 Stationary marked point processes 

An MPP Xm = Y2 n >l 3{X n ,M n ) * s a ran dom locally finite counting measure acting on the Borel 
sets of R d X M with atoms (X n , M n ) , where the marks M n belong to some Polish mark space 
M endowed with the Borel cr-algebra B(M). Throughout we assume that Xm is simple, i.e. 
all locations X n in R" have multiplicity 1 regardless which mark they have. Mathematically 
spoken, Xm is a measurable mapping Xm ■ ^ — > Nm from some probability space (f2, ^4, P) 
into the set Nm of counting measures </?(•) on B(R d x M) satisfying tp(B x M) < oo for 
all bounded B £ Z3(R d ) , where Nm is endowed with the smallest cr-algebra Mm containing 
all sets of the form {ip 6 Nm : ip(B X C) = j} for j > 0, bounded B G B(R d ) , and 
C € B(M) . In what follows we only consider stationary MPPs, which means that the 
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distribution Px M (') = ^{-^m £ (*)) °f -^M 011 -A/m is invariant under location shifts of the 
atoms, i.e., 

X M = S (X n -x,Mn) for a11 x£R d . 
n>l 

Provided that Xm is stationary and the intensity A = EX^-QO, l) d x M) is finite we have 
EX M (B x C) = \B\ EX M ([0, l) d x C) for all bounded B G B(R d ) and C G B(M) , where | • | 
denotes (i-dimensional Lebesgue measure. 

2.1 Palm mark distribution 

For a stationary MPP the probability measure P M 011 B(M) defined by 

P M (C) = jEX M ([0, l) d xC), C G S(M) , (2.1) 

is called the Palm mark distribution of Xjvf. It can be interpreted as the conditional distri- 
bution of the mark of an atom of Xm located at the origin o . A random element Mq in M 
with distribution P M is called typical mark of Xm- 

Definition 2.1. An increasing sequence {W&} of convex and compact sets in W 1 such that 
Q(Wk) = sup{r > : B(x,r) C W k for some x G W k } — > 00 as k — > 00 is called a convex 
averaging sequence (briefly CAS). Here B(x,r) denotes the closed ball (w.r.t. the Euclidean 
norm || • ||) with midpoint at x G M. d and radius r > . 

Some results from convex geometry applied to CAS {Wfc} yield the following inequalities 

_ L_ < H*-i(dw k ) < _J_ and \w k n(w k -x)\ < dM 

g(W k ) ~ \W k \ ~ g(W k ) \W k \ - g(W k ) K ' 

for < g(W k ) , where T-Ld-i(dW k ) is the surface content (i.e. (d— l)-dimensional Hausdorff 
measure) of the boundary dW k , see [3] and [16] for details. 

If Xm is ergodic (for a precise definition see [5], p. 194), the individual ergodic theorem 
applied to MPPs (see Theorem 12. 2. IV and Corollary 12. 2. V in [5]) provides the P — a.s. 
limits 

for any C G B(M) and an arbitrary CAS {W k } . 



2.2 Factorial moment measures and the covariance measure 

For any integer m > 1, the mth factorial moment measure ct^ of the MPP Xm is defined 



in 



on B({R d x M) m ) by 

a£(x(5 l xC ! ))=E ^ ]{l Bi (X ni )l Ci (M ni )) , (2.4) 
i=i '—^ 
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where the sum n m >i runs over an m-tuples of pairwise distinct indices m, . . . , n m > 1 

for bounded Bi G B(R d ) and Cj G £>(M) , i = 1, . . . , m. We also need the mth factorial 
moment measure ct^ of the unmarked PP X(-) = Xm((-) X M) = ^ n >i ^x n (') defined on 
B{(R d ) m ) by 

a { x ] (x^B-\ = a { ™l ( x (B t x M) ) for bounded B u . . . , B m G £(M d ) . 

The stationarity of Xm implies that a x is invariant under diagonal shifts, which allows to 

define the mth reduced factorial moment measure otxted uinc l uer y determined by the following 
desintegration formula 

4° (.x B,)=\f ag ed ( x ( Bl -x))dx. (2.5) 
v=i / Jb 1 »=2 

We need a condition of weak dependence between parts of X defined over distant Borel sets 
which can be expressed by the (factorial) covariance measure on B{(R d ) 2 ) defined by 

7^(51 x B 2 ) = a ( ${B 1 x B 2 ) - A 2 \B 2 \ . 
The reduced covariance measure 7^ ■ : S(M d ) — )• [—00,00] is in general a signed measure 

(2) (2) 

defined by (2.5) with ^ x instead of a x , which means that 

lx]red( B ) = a x\ed( B ) ~ X \ B \ for bounded B G B{R d ) . 

For more details on factorial moment measures and measures related with them we refer to 
Chapters 8 and 12 in [5]. 



2.3 ra-point Palm mark distribution 

For fixed mark sets C\, . . . , C m G £>(M) , m > 1 , the mth factorial moment measure 
defined by (2.4) can be regarded as a measure on the Borel sets B((R d ) m ) depending on 
C±, . . . ,C m . This new measure is absolutely continuous w.r.t. the mth factorial moment 
measure cS x \ Thus, the Radon-Nikodym theorem (cf. [7], p. 90) implies the existence 

of a density Pf}''"' Xm (Ci x ••• x C rn ), which is uniquely determined for -almost all 
(xi, . . . , x m ) G (R d ) m , such that for any B±, ... , B m G B(R d ), 

a™(x(B i xC i )) = ! P x M ^ Xm ( x a) a i x a) (d(x 1 ,...,x m )). (2.6) 

i=i 

Since the mark space M is Polish, this Radon-Nikodym density can be extended to a regular 
conditional distribution of the mark vector (Mi, . . . , M m ) given that the corresponding atoms 
Xi, . . . , X m are located at pairwise distinct points x±, . . . , x m , i.e., 



P^-' Xm (C)=¥((M 1 ,...,M m ) eC\X 1 = x 1 ,...,X m = x m ) for C G £>(M m ) . 

This means that the mapping (x±, . . . , x m , C) 1— > Pf^''"' m (C) is a stochastic kernel , i.e., 
P^'-' Xm (C) is £((R d ) m )-measurable in (x 1 ,...,x m ) G {R d ) m for fixed C G B{M m ) and a 
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probability measure in C G B(M m ) for fixed (x\, . . . ,x m ) G (JR^)" 1 . For details we refer to 
[18], p. 164. The regular conditional distribution P^'-' Xm (C) for C G B(M m ) is called the 
m-point Palm mark distribution of Xm ■ This stochastic kernel is only of interest for m-tuples 
{x\, . . . , x m ) of pairwise distinct points Xi G R d , i = 1, m. In case of a stationary simple 
MPP Xm it can be shown that 



M 



(C) = p^-*i...,x m -zi( C ,) for c G ^(M m ),m>l 



and any xi, . . . , x m G R d with x, 7^ Xj for i ^ j ■ In this way the Palm mark distribution 
defined in (2.1) can be considered as one-point Palm mark distribution. 

The following result is crucial to prove asymptotic properties of variances estimators of the 
empirical mark distribution. It generalizes an analogous result stated for unmarked PPs in 
[17] to MPPs by involving the notion m-point Palm mark distribution for m = 2,3,4. The 
proof is just a slight extension of the one of Lemma 5 in [17] by using the relation (2.6) for 
m = 2, 3, 4 . The details are left to the reader. 

Lemma 2.1. Let Xm = Y^ n >i ^(X„,M„) be an MPP satisfying E(Xm(B x M)) 4 < 00 for all 
bounded B G B(R d ), and let f : R d x M. d x M 2 i-> R be a B or el-measurable function such that 



the second moment of q>l \ f(X p , X q , M p , M t 



exists. Then, 



Vax^jf f(X p ,X q ,M p ,M q ) 
p,<?>i 



(2.7) 



f{xi,x 2 , ui, u 2 ) f{xi,X2,u 1 ,u 2 )+f{x 2 , xi,u 2 , ui) Pm' X2 (d(«i, «2))ax {d(x 1 ,x 2 )) 



+ 



J J f(xi,x 2 ,ux,u 2 ) f(x 1 ,x 3 ,ui,u 3 ) + f(x 3 ,x 1 ,u 3 ,u 1 ) 



+ 



+ f(x2, x 3 ,u 2 ,u 3 ) + f{x 3 ,x 2 ,u 3 ,u 2 ) P^' X2 ' X3 (d(ui,u 2 ,u 3 ))a ( -x > {d(xi,x 2 ,x 3 )) 



(R d ) 4 M 



j /(*!, 

4 M 4 



X 2 ,Ul,U2)f(x 3 ,X4,U 3 ,U4,] 



M 



(d(ni, u 2 , u 3) ii4))a^ (d(x 1 ,x 2 ,x 3 ,x i )) 



P%' x * (d(m, u 2 )) P x t x * (d(« 3 , «4)) a? (d(x 1? x 2 )) a£> (d(x 3 , x 4 )) 



.,(2) , 



2.4 /3-mixing coefficient and covariance inequality 



For any B G i3(R ), let Ax M (B) denote the sub-c-algebra of A generated by the restriction 
of the MPP Xm to the set BxM. For any B, B' G B(M. d ) a natural measure of dependence 
between Ax M (B) and Ax M (B') can be formulated in terms of the (3— mixing (or absolute 
regularity, respectively weak Bernoulli) coefficient 



p{a Xm (b),a Xm (b')) = \ sup y,\ n 4) - p ( a op(40 



(2.8) 



where the supremum is taken over all finite partitions {A{\ and {A'j} of 0, such that G 
Ax M (B) and G .Ax M (-£>') for all i,j , see e.g. [6], [12] or [26]. It should be noticed that the 
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supremum in (2.8) does not change if the sets and A 1 - belong to semi-algebras generating 
Ax M (B) and Ax M {B'), respectively. To express the degree of dependence of the MPP Xm 



for disjoint sets K a = 
functions /3^ M , P x 



-a, a] d and K£, b = M. d \K a +b, where b > 0, we consider non-increasing 



i,oo) 



0, oo) such that 



for i < a < b , 



1 for \<b<a 



(2.9) 



A stationary MPP Xm is called /3-mixing or absolutely regular, respectively weak Bernoulli if 
both (3-mixing rates (3 Xu (r) an( ^ @x ( r ) tend to as r — > oo. By standard measure-theoretic 
approximation arguments it is easily seen that any stationary /3-mixing MPP Xm is mixing 
in the usual sense and therefore also ergodic, see Lemma 12. 3. II and Proposition 12. 3. Ill 
in [5] Vol. II p. 206. In order to prove CLTs we need further conditions on the decay of 
the /3-mixing rates (3* x (r) and (3* Xhi (r) on the right-hand side (rhs) of (2.9). For this we 
formulate 

Condition (3(d): There exists some 5 > such that E(X M ( [0, l] d x M)) 2+S < oo , 
r d - l (P XM (r)) &,{2+&) dr < oo and r 2d ~ 1 (3 Xu (r) — » . 

The following type of covariance bound in terms of the /3-mixing coefficient (2.8) was first 
stated in [26]. 

Lemma 2.2. Let Y and Y' denote the restrictions of theJAPP X M to B x M and B' x M 
for some B, B' £ B{M. d ) , respectively. Furthermore, let Y and Y' be independent copies of 
Y and Y' , respectively. Then, for any Mm ® A4«-measura6te function f : Nm x Nm — > [0,oo) 
and for any i] > 

\Ef{Y,Y')-Ef{Y,Y')\ < 2(3(A Xm (B),A Xm (B'))^ 

x maxjtE/^^Y')) 1 ^ (E/ 1+ ^(Y, Y')) ^ } • (2-10) 

If / is bounded, then (2.10) remains valid for r/ = oo . In the particular case f(y,y') = 
fi(y) f2(y') an d r\ = 5/2 for 5 > 0, the Cauchy-Schwarz inequality applied to the expectations 
on the rhs of (2.10) yields 

| Cov( MY), h{Y') ) | < 2 || h{Y) \\ 2+5 || h{Y') \\ 2+s {(3(Ax m (B),A Xm (B') , ( 2 .11) 
where || Z \\ q = (E|Z|«)V<z is the L« —norm (q > 1) of a random variable Z . 

3 Results 

3.1 Central limit theorem 

We consider a sequence of set-indexed empirical processes { Yfc(C) , C £ £>(M) } defined by 

y fc (C) = -±= £ iw k (Xn){i c (M n )-p M (c)) = {(P M ) k {C) - P° M (C)) , (3.1) 
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where {W^} is a CAS of observation windows in M d . We will first state a multivariate 

CLT for the joint distribution of Y k {C\), . . . , Y k (C(). For this, let t —^' denote convergence 
in distribution and Mi{a, X) be an ^-dimensional Gaussian vector with expectation vector 
a£l ; and covariance matrix I] = (<Xjj)| =1 . 

Theorem 3.1. Let Xm be a stationary MPP with A > satisfying Condition /3(5). Then 

Y k = (y fc (d),...,y fc (C,)) T A JVi(o<,E) /orany Ci,...,QeB(M), (3.2) 
where oi = (0, . . . ,0) T and i/ie asymptotic covariance matrix S = (<7ij)| - =1 is given fry t/te 

ay = lim EnCCOnCC,-). (3.3) 

The above result can also be stated in terms of the empirical set-indexed process {Z k (C), C 6 
B(M)}, where 

£*(C7) = ( Afe ) _1 Yfc(C) = vW((^) fc (C) - P? f (C)) . 

In other words, as refinement of the ergodic theorem (2.3), we derive asymptotic normality of 
a suitably scaled deviation of the ratio-unbiased empirical Palm mark probabilities (-PJu-)fc(C) 
from P M (C) defined by (2.1) for any C £ £>(M) . Since Condition f3(5) ensures the ergodicity 
of Xm , the first limiting relation in (2.3) combined with Slutsky's lemma yields the following 
result as a corollary of Theorem 3.1. 

Corollary 3.2. The conditions of Theorem 3.1 imply the CLT 

Z k = (Z k (Ci), Z k (C e )) T A M e (o e , A" 2 E) . 

k— >oo 



3.2 /3-mixing and integrability conditions 

In this subsection we give a condition in terms of the mixing rate /3 Xm (r) which implies finite 

(2) 

total variation of the reduced covariance measure ^) x ref ^ and a certain integrability condition 
(3.5) which expresses weak dependence between any two marks located at far distant sites. 
Both of these conditions are needed to get the asymptotic unbiasedness resp. Inconsistency 
of some estimators for the asymptotic covariances (3.3). 

(2) (2) 

Note that the total variation measure |7^- re( J of 7 X red * s defined as sum of the positive part 

(2)+ (2)— (2) 

iX red ano - negative P ar t 7x red °f ^ ne J° rc ian decomposition of "fx red> i-e., 

J 2 ) _^(2)+ _^(2)" and | (2) | (2)+ (2)- 

IX, red ~ IX, red 'X,red dnu I 'X,red\ ~ 'X.red ' 'X,red > 

(2)-\- (2) 

where the positive measures "f x red ano - iX red are mu tually singular, see [7], p. 87. 
Lemma 3.1. Let Xm be a stationary MPP satisfying 

POO 

E(X M ([0,l] d x M)) 2+5 < oo and / r^ 1 ((3 XM (r)) 5/{2+5) dr < oo for some 5>0. 



(2) ,1 

Then "fx red ^ as fi n tt e total variation on K , i.e., 

(2) 
fx,red 



l7?Ll(M d )<oo. (3.4) 



Furthermore, for any C\^C 2 E S(M) 

/ Ptf(C x x C 2 ) - i^(Ci) P M (C 2 ) a%] red {dx) < oo . (3.5) 
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3.3 Representation of the asymptotic covariance matrix 

In Theorem 3.1 we stated conditions for asymptotic normality of the random vector Y&. 
Clearly, (2.1) and (3.1) immediately imply that EY k (C) = for any C £ B{M). The following 
theorem gives a representation formula for the asymptotic covariance matrix 13 . 

Theorem 3.3. Let Xm be a stationary MPP satisfying (3.5) and let {Wk} be a CAS. Then, 
the limits in (3.3) exist and take the form 

*ij = A(P A / (C i nL7 i )-P A / (C7 i )P A / (C,))+A / (PtfiGxCj) (3.6) 

- p*r(G x m) pu^) - Ftf{Cj x m) p^a) + mci) mcj) ) 4]^) ■ 

In particular, if Xm is an independently MPP, then 

aij = A( P M (d n Cj) - P°f(Ci) PSi(Cj) ) . (3.7) 



3.4 Estimation of the asymptotic covariance matrix 

In Section 6 we will exploit the normal convergence (3.2) for statistical inference of the typi- 
cal mark distribution. More precisely, assuming that the asymptotic covariance matrix 53 is 
invertible, we consider asymptotic x 2 -g°°dness-of-fit tests, which are based on the distribu- 
tional limit 

Y k Y k , > Xli 

which is an immediate consequence of (3.2) and Slutsky's lemma, given that is a consistent 
estimator for S. Here we use the notation Y k = (ifc(Ci), . . . , Yfc(C^)) T (see (3.1)) and the 
random variable xj is X 2 ~distributed with I degrees of freedom. In the following we will 
discuss several estimators for 53. Our first observation is that the simple plug-in estimator 
= (Yjfc(Cj)Yfc(Cj))^ ._ 1 for S is useless, since the determinant of xjj. vanishes. Instead 

of 53^ we take the edge- corrected estimator = ((o^ 1 =1 with 



(a£\ = ^Y.^w k {X p )(K Crr[C] (M p )-PUCi)P° M {C 3 )) (3.8) 
P>1 

^ l Wk (X p )l Wk (X q )(l Ct (M p ) - P° (Q) ) {l Cj (M q ) - P° M {C 3 )) 



+ E 



\(w k -x p )n(w k -x q 



As an alternative, which can be implemented in a more efficient way, we neglect the edge 

T {2) ) Y 



correction and consider the naive estimator = =1 f° r 53 with 



= r ^E I ^( x p)( I ^( M p)- p A/(a)p,°/(c J )) 

+ iiE^ l W /J^)l^(A g )(l Cl (M p )-P^(C i ))(lc J (M g )-P A ° / (C J )). 

p,q>l 
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Theorem 3.4. Let Xm be a stationary MPP satisfying (3.5) and let \Wk\ be a CAS. Then 
ffiij Jk ^ an un biased estimator, whereas (3y is an asymptotically unbiased estimator for 

i h 3 = 

Remark: In general, neither )fe nor (3^ ) fe are L 2 -consistent estimators for cr^, even if 
stronger moment and mixing conditions are supposed. 

According to Lemma 3.1, the integrability condition (3.5) in Theorems 3.3 and 3.4 can be 
replaced by the stronger Condition /3(S). In order to obtain an L 2 -consistent estimator, we 
introduce a smoothed version of the unbiased estimator in (3.8), which is based on some 
kernel function and a sequence of bandwidths depending on the CAS {Wfc}. 

Condition (wb): Let w : R h- > K be a non-negative, symmetric, Borel-measurable kernel 
function satisfying w(x) — > w(0) = 1 as x — > 0. In addition, assume that w(-) is bounded 
by m w < oo and vanishes outside B(o,r w ) for some r w G (0, oo). Further, associated with 
w(-) and some given CAS {Wfc}, let {b k } be a sequence of positive bandwidths such that 

>b k — ►() , b d k \W k \ -^oc and fcf d |W fc | — > 0. (3.9) 



2dr w \W k \ 1 / d ~ fc^oo ' fc 1 'fe^oo k 1 'fe- 



Theorem 3.5. Let {W k } be an arbitrary CAS andw(-) be a kernel function with an associated 
sequence of bandwidths {b k } satisfying Condition (wb). If the MPP Xm satisfies 



E(A M ([0,l] d xM)) 4+5 <oo and f°° r^ 1 (P* Xm (r)) 5/(4+5) dr < oo (3.10) 
for some 5 > , then 



, ('*-*'h i », 
where (&■>■ ) k is a smoothed covariance estimator defined by 

{a^X = -^^l Wk (X p )(l c% ncAK)-PM(a)PMCj)) 



+ !^ W !^ TO ( ^ W ~ p m( c -) ) ( W ~ ) ^yii^g- x p 



|(w fe -x p )n(w fe -x 9 )| V 6 fc |w fc |Vd 

Remark: The full strength of condition (3.10) on the /3-mixing rate (3 X (r) introduced in 
(2.9) is only necessary to prove the consistency result of the preceding Theorem 3.5. However, 
the /3-mixing rate /3 Xm (t) m Condition (3(5), which is needed to prove (3.4) and (3.5) as well 
as Theorem 3.1, can be defined by the slightly smaller non-increasing /3-mixing rate function 

/3* XM (r)=f3(Ax M (Ka),Ax M (K c a+ r)) for r>o=l/2. (3.11) 

Moreover, in order to prove Theorem 3.1, condition (3(5) relying on the /3-mixing coefficient 
considered in (2.8) with j3 x (r) and Px M (r) given in (3.11) and (2.9), respectively, can be 
relaxed by using the slightly smaller a-mixing coefficient 

a(A XM (B),A XM (B')) = sup{ \F(A n A') - F(A)F(A')\ : A G A Xm (B),A> G A Xm (B')} 
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instead of (2.8). The corresponding a-mixing rates a* XM {r) and ot*^ M {r) are then defined 
in analogy to (3.11) and (2.9), respectively. A covariance inequality for the a-mixing case 
similar to (2.11) can be found in [6], see [14] for an improved version. Despite of the subtle 
differences between the discussed mixing conditions, we prefer to present our results under 
the unified assumptions of Condition /3(5) and (3.10) with fi^ (r) as defined in (2.9). It 
seems to be difficult to identify models where these differences are relevant. 



4 Proofs 



4.1 Proof of Theorem 3.1 

By the Cramer- Wold technique, the multivariate CLT stated in (3.2) is equivalent to 

s T Y k = siY fe (Ci) + ... + s t Y k [C t ) A A/i(0, a 2 ) with a 2 = s T 'Es (4.1) 

k— >oo 

for any s = (s 1; . . . , sg) T £ R f / Of , 

To prove (4.1) we extend Bernstein's classical blocking method for weakly dependent random 
fields over a cubic index set of Z d , see e.g. [4], [12] or [21], to /3-mixing fields indexed by 
elements of H k = {z G 7L d : E z C W k } , where E z = [-1/2, l/2) d + z for z G 7L d and {W k } is 
an arbitrary CAS. The proof of (4.1) is divided into four steps. 

Step 1. Bounds and asymptotics for the variance of the sum 

In view of (3.1) we may write 

s^Y k = ^={V k + Vl), where V k = £ U z and V^= £ 

with 

TJW =J2 I E*nW k (Xn)9(M n ) , U z = Y,^E z {Xn)g{M n ) for z G Z d , 

n>l n>l 

3H k = {z G 7L d \ H k : | E z n W k \ > } and g(M n ) = £- =1 Si(l Ci (M n ) - P^Q)) . Clearly, 

EU { z k) = EU Z = and max{ | U { z k) | , | U z | } < c(s)X M (£* x M) for z £ Z d , since |s(M n )| < 
c(s) = |si| + • • • + |s^| . Hence, by stationarity of Xm and Condition /3(5), 



Lax{|| C/ Z (fc) \\ 2+S , || E7, ||2+d} < c(s) ||A M ([0, l) d x M)|| 2+^ for z G n ' 



In the following we use the maximum norm \z\ = maxi<j<^|zi| to express the distance of 
z = (zi,...,Zd) G Z rf to the origin o. By applying the covariance inequality (2.11) together 
with Condition (3(8) , we obtain 

Var(V^) = E U v k) U z k) ^ Yl E \ U i k) U ^ I - # dHk Yl E ^ U ^ I 

y,z&dH k y,zedH k z ez d 



< 2c(s) 2 \\X M ([0,l) d x M)||2 +(5 #aff fc ^ (/3(^x M (£oMx M (£ 2 )) 



2+« 
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<2c( S ) 2 ||XM([0,l) d xM)||2 +5 #^(3 d + £ fe(N-l))^ 

z: \z\>2 

< 2c(s) 2 ||X M ([0,l) d x M)\\ 2 2+s #dH k (3 d + 2d ^(2n + 3) rf " 1 C 4 ' 2 ) 

n>l 

for some constant c\ = c\(s, d,5) > where the relation #{z £ Z rf : |z| = n} = (2n + l) d — 
(2n — l)" 1 < 2 d (2n+ has been used. A simple geometric argument shows that each unit 

cube E z hitting the boundary dW k is contained in the annulus dW k © B(o, yd) implying 
that 

#dH k < \dw k © b(o, Vd)\<2( \w k © b(o, Vd)\ - \w k \ ) . 

Steiner's formula (cf. [24], p. 600) applied to the convex body W k reveals that the volume 
\W k © B(o, yfd)\ — \W k \ does not decrease when W k is replaced by a larger convex body, 
e.g. by d 3 / 2 R k from relation (4.8) below, where the hyper-rectangle R k has edge lengths 
. . . , ajj*^ . Replacing additionally B(o, \fd) by the cube [—y/d, \/d\ d we get 

\W k ®B(o,Vd)\-\W k \ < \d^ 2 R k ®[-Vd,Vd] d \-\d 3 / 2 R k \ 



d 



i=l 



2 s^y. «?>■■■ «e?.Wi+ 3) ■■■(«?'+ 2 

2 



d/ V « d 



< 2 d - 1 d( M - 2 )/ 2 ^_ 1 (5 J R fc ), if minaf^ 

l<i<d a 

Hence, since (4.8) implies U d _i(dR k ) < Hd-i(dW k ) and d 3 / 2 mini<j<^ a$ > 2 g(W k ) , it 
follows that #dH k < 2 d d^ d ~ 2 ^ 2 H d -i(dW k ) if g(W k ) > y/d, which in turn by combining 
(2.2), (4.2) and the inclusion #H k < \W k \ < #H k + #dH k implies that 

Y^m iC2 n^m^ and (43) 

I I I Wit I e(Wik) Wifel fc->oo 



for any CAS {W^} . Thus, by a standard Slutsky argument, (4.1) is equivalent to 

AMOV). (4.4) 



V k D, 



The technique used above to estimate Var(V^) will in the following be applied to show that 
a 2 = lim Var( S T r fc ) = lim = TE(U U Z ). (4.5) 

k— too k— ¥00 #tlu 

W K Z&i 

The series on the rhs of (4.5) converges absolutely as immediate consequence of the estimate 
Var(14) < #H k \E(U U Z )\ < Cl #H k , 
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where the positive constant c\ is the same as in (4.2). The Cauchy-Schwarz inequality and 
the previous estimates of Var(Vfc) and Var(V fc ') show that 



Var(s 1 Y k ) 



Var(y fc ) 



< 2 \CoHVk,V0\ Var(^) 2c^#H k #dH k q#dH k 
I W fc | |W fc | - " |Wfc| " |W fc | 



proving the second equality in (4.5). To prove the third equality in (4.5) we use the identity 



1- £ E( W= S *< g '2*'-'» E (I ^. 



and the geometric inequality (following from the very definition of H k and dH k ) 

#(H k n(H k - z ))<\w k n(w k -z)\< #(H k n(H k - z )) + #dH k + #d(H k - z) 

for z £ Z rf . This fact combined with (2.2) and (4.3) shows that 

#(H k n (gfc - z)) d 

— — > 1 for any fixed z£Z 

^fH k k— ¥oo 

proving the third equality in (4.5) by applying the dominated convergence theorem. 
Step 2. Passage to bounded random variables by truncation 

For any fixed a > we define the random field {U z (a) , z £ H k } of the truncated (and 
centered) random variables and the sum V k (a) by 

U z (a) = U z l {lUzl < a} -E(U z l mi < a} ) and V k (a) = ^ U z (a) (4.6) 

z£H k 

so that, for any z € H k , 

\U z (a)\<2a and ( E| U z - U z (a) \ 2+s ) = \\ U a - U Q (a) || 2+<5 — > . 

a— >oo 

By quite the same arguments as used in Step 1 based on the covariance inequality (2.11) and 
Condition (3{5) , we find that 

Var( V k -V k (a))<2 #H k \\ U Q - U (a) f 2+s ( 3 d + 2 d ^(2?i + 3)^ (/% M (n)) ^ ) 

n>l 

for k > 1 . Hence, by Slutsky's lemma, the weak limits of V k /\/#H k and V k (a)/\/#H k as 
A; — 7- oo are arbitrarily close whenever a > is large enough. It therefore remains to prove 
the CLT in (4.4) for the bounded random variables in (4.6), i.e., for any fixed a > , 

iM^L A A/i(0, ct 2 (o)) with ( r 2 (a)= VE(/ (a)[/ 2 (a). (4.7) 
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Step 3. Adaptation of Bernstein's blocking method to non-cubic index sets 



We start with some preliminary considerations. A well-known result from convex geometry 
first proved by F. John, see e.g. [1], asserts that there exists a unique ellipsoid £k (called 
John ellipsoid) of maximal volume contained in Wk with midpoint c(£ k ) and semi-axes of 
lengths , . . . , such that £ k C Wk Cj c(£ k ) + d(£k — c(£ k ) ) ■ 

Further, it is easy to determine a unique hyper- rectangle Rk centered at the origin o cir- 
cumscribed by £k — c(£k) with edge-lengths = 2e[ k ^/y/d for i = l,...,d such that 
£k — c(£k) C y/dRk and finally 

R k Q£ k - c(£ k ) C W k - c{£ k ) Cd(£ k - c(£ k ) ) C d^ 2 R k . (4.8) 

Since the MPP Xm observed in the CAS {Wjfc} is stationary, we may assume that c{£ k ) = o 
and without loss of generality let the edge lengths of Rk be arranged in ascending order 

(k) (k) 

a\ < • • • < a\ (possibly after renumbering of the edges). Note that Rk is not necessarily 
in a position parallel to the coordinate axes. But there is an orthogonal matrix Ok such that 

(4.9) 



* r «<*> a[ k) 



OkRk — x 
i=l 

Let {pk} and {qk} be two sequences of positive integers (which will be specified later) satis- 
fying pk > qk — > oo and qk/pk — > 0. We define two types of pairwise disjoint cubes 

k— >oo k— >oo 

pW =Pi k) + (2p k + q k + l)y and Q<*) = Q (k) + (2 Pk + q k + 1) y for y G Z d , 

where = {-p fe , . . . , 0, . . . ,p k } d and = {-p k , . . . , 0, . . . ,p k + g fc } d for fe > 1 . 

Now, we describe how to modify Bernstein's blocking method in order to prove the CLT stated 
in (4.7). For the family of block sums 



v i k) ( a )= E u *( a ) for y£G k = {zez d :P^nH k ^®} 

we shall show in Step 4 that 

by assuming the mutual independence of the random variables (a) , y £ G k , which can 
be justified by Condition /3(5). Moreover, it will be sufficient to prove (4.10), since we will 
show below that 

-^=(v k (a) - ]T f(a)) A 0. (4.11) 

Next, we specify the choice of pk and qk in dependence on the edge lengths of Rk and 
the supposed decaying rate of Px M ( r ) ■ ^ n y i ew °f Q(W k ) — > oo and (4.8) it follows that 
min{a^, . . . , a^} — > oo as k — > oo . Note that the choice p k = [ £ fc | Vl/fc 1 1 / 2d J as in case 
of a cubic observation window with a certain null sequence {£&} , see [12], does not always 
imply (4.10) and (4.11) if at least one of the first d — 1 ordered edge lengths of Rk increases 
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(k) 

very slowly to infinity So one has to choose pk large enough but much smaller than a K d ' . For 
this purpose put ffc(s) = ( • . . . • ) 1 /( 2d ~ s ) for each s G {0, 1, . . . , d — 1} . Because of 
1 fix ( r ) — ^ as r — )■ oo, there exist non- increasing sequences £fc(s) of positive numbers 



such that 



2d-l 



£k(s) — > 0, £fc(s)r fc (s) — ^ oo, and — P*x M ( £ k(s)r k (s)) — >0. (4.12) 

k—>oo k— >oo £k\ s ) k— >oo 



Let £fc = max{efc(0), . . . , £k{d— 1)} and Pk{s) = ^ fk( s ) an d select Sk to be the smallest 

(k) 

number s G {0, 1, d — 1} such that a),^ > 2pk(s) + 1 for k > ko , where /co is a sufficiently 
large positive integer. Thus, we define the integer sequences pk and qk by 



Pk = [Pk{sk)\ = L e fc Sfc) ^fc(sfc)J and % = Lefcr fc (s fc )J for k > k . 



(4.13) 



Further, we need lower and upper bounds for the number iVfc of cubes Qy = [—pk — \-,Pk + 

% + i) 01 + (2p k + q k + l)y hitting F fe , i.e., N k = #{y G Z d : n H k ± 0} . For this put 

N^icw) = #{y G Z d : Q { y k) n (Lf\c) + w) ^ $} for io G M d and some real c> 0, where 

Lf\c) = Ol{(xi,...,x d ) G R d : -caf ] /2 < Xj < caf ] /2, x< = for i / j} . The 

following rough estimates of iVj (c, u;) from below and above can be obtained by elementary 
geometric arguments: 



caf/y/d 
2pk + Qk + 1 



+ 1 < N^ k \c,w) < d 



ca 



(fe) 



2pfe + % + U 



+ 2 I for any u; G 



Hence, by (4.8) and (4.9) the minimal number N^J and the maximal number N^ax of cubes 
Q-y hitting ff^ satisfy the inequality 



n(l 



aa\ L2p fc + g fc + l 



+ l) < iV« < iV fc < tf<W < 



2pfc +qk + 1 



+ 2 



In view of the above choice of s = Sk and (4.13), the number Nf. allows the estimate 

,( fc ) .J k ) J k ) J k ) 



C3 



° Sfe+1 " ' ° d < iVfe < c 4 Qsfc+ V ' ' Qd for all k > k 



d—Sh 
P k 



(4.14) 



Pk 



with positive constants 03,04 only depending on the dimension d. Combining the obvious 
fact that #G k < N k with (4.12), (4.13) and (4.14) (with p k > 1 and e k < 1) we arrive at 



d— 1 a** ^ ^ ( r fc( s fc)) o** / \ n 

■ Px M yi k > , — > u • 



l/2d 



k~ ¥00 



Likewise, by (4.8) and < 2p k + 3 for i = 1, s k , 



d—Su s k 

Pk TT (k) , 

v 7 ^ f = \ 3 ~ 



Pk 



rk(sk) 



(2d-s k )/2 



< c 5 y/ek~ — > . 

k— >oo 



(4.15) 
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Finally, we show that 



l — E #(f nff t ) — >i, (4.16) 



which, by the results of Step 1, is equivalent to 

±- ElK'Vf 1 )""'.!^, (4.i7) 



where Py fe ^ = [—p k — \,Pk + \) d + (2pfc + Ofc + 1) y for y £ Z, d . To estimate the volume of 

~(k) (k) r l 

the space in W k outside the union of cubes Py we introduce equidistant slices in M of 
thickness and distance 2p k + 1 defined by 

S\f = {( yi ,...,y d )eR d : (2j + 1) (p k + 1) + j q k < Vi < (2j + 1) (p k + I) + (j + 1) 9fc } 

for i = 1, ...,d and j £ Z 1 . By (4.8), (4.9) and the choice of and q k it might happen that, 
for at most s k coordinates i £ {1, d}, 5^ D VFfc = for all intergers j . For the remaining 
coordinates % G {1, d} there exist sequences of integers n k {i) (at least one of them tends to 
infinity as k — > oo) such that D -Rfc ^ for |j| < nfc(i) (and 5^ D = for |j| > n k (i)) 
and 

E I S if n d3/2 ^fc I < c 6 - for k > k , 
I ^fe I I v Pk 

where C6 depends only on d. This estimate and the evident inequalities | R k \ < \ W k | and 

E I ( Ql k) \ ?i k) ) n ^l<E E 1 n d3/2 I 

show that the lhs of (4.17) is bounded by a constant multiple of q k /p k so that (4.16) is finally 
proved by (4.13). 

Step 4. Approximation by sums of independent random variables 

For brevity put P k = UyeGfcC^i/ H -£/&) . Again by applying the covariance inequality (2.11) 

and Condition f3{5) to the stationary random field {U z (a),z € H k } (with |C/ z (a)| < 2a and 
thus <5 = oo), we find in analogy to (4.2) that 

jL-E(v H (a)- E^W) 2 = IF" E E(C/»£Ua)) 

< 8 a 2 ( 3 d + 2 d E( 2 " + ^ W ) * {H ^ ^ • 

n>l * k 

From (4.16) it is immediately clear that the ratio in the latter line disappears as k — > oo, 
which confirms (4.11). Thus, in view of Slutsky's lemma, it remains to prove (4.10). We will 
do this under the assumption of mutual independence of the block sums V y (a) , y G G k . For 
this reason we show that the characteristic function Eexpjff ^2 y&Gk ^2/ ( a )} differs from the 

product Hj/eGfe Eexpjit V y (a)} uniformly in t E I 1 by certain sequences tending to zero as 
k — > oo . 
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Setting n k = #G k we may write 

Cj = exp{r£ V^\a)} for yj £ G k with j = 1, . . . , n k 
Using the algebraic identity 



n fe n k 



nt— 1 



E II Hi II E ^ = E E ^ • • ■ E ^-l ( E ^i+l ■ ■ ■ in k - Efe +1 • • • £n fc ) ) 



3=1 3=1 3=1 

and < 1 for j = 1, we get 



ra/t-i 



Eexp{it £ F«(a)}- J] Eexp{ityW(a)}| < £ | Cov(^, & +1 • • • £nj 
?/GG fe yeG fc i=l 



By the stationarity of Xm we may assume that the real as well as the imaginary part of £j is 
measurable w.r.t. the a— algebra A.x M {Kp k +i/2) and the product ■ ■ ■ £ nfe is measurable 
w.r.t. Ax M {^. d \ K pk+qk+ i/ 2 )- By applying the covariance inequality (2.11) with 5 = oo (to 
the real and imaginary part of resp. ■ ■ ■ £,n k ) and using (2.9) we find that 

Covfo^+i-.-eol < 8P{Ax M (K Pk+ i/2),Ax M (^ d \K Pk+qk+ i/2)) 

< 8 + 1/2)^ 



Since = #Gfe < Nk it follows with (4.14) that 



sup 



Eexp{it V v (k) ( a )}- II Eexp{ityW(a)} < 8n fc (p fc + 1/2 ^ ^ M fe) 



k— >oo 



0. 



The latter relation and the Berry-Esseen bound in the CLT for independent random variables 
(which can be expressed by the third-order Lyapunov ratio, see e.g. [4], p. 204, and references 
therein) reveal that (4.10) holds if 



4 fc) (a) 



(^(«)) 3/2 



£E|vf)(a) 



yeG k 



| 3 — > and 

k— >oo 



a 2 (a), (4.18) 



where c 2 (a) is defined by (4.7) and c 2 (a) = X^eGfc E( (a) ) 2 coincides with the variance 
of T4(a) in case of independent block sums Vy (o), y 6 . 



It is easily seen that | (o) | < 2 a # (P y W HH k ) <2a (2p k + l) d "^ Tli=i( rf3/2 a f J + x ) 



>(*) 



(fc) 



and therefore 



Lf(a)<2an(d 3/ \ (fc) + l) 



(fe) 



(2p fc + 1 



\d—si 



o „d—s k s k 

2a Pk TT„( fc ) 



i=l 



£^M>*<c 7 ^II«! 



yeG fc 



i=i 



with some positive constant C7 only depending on d. In combination with (4.15) the second 
relation in (4.18) for cr 2 (a) > yields Lg (a) — )• 0. Hence, the first part of (4.18) is proved. 

k— ¥00 

To accomplish the proof of (4.18) we remember that cr 2 (a) is the asymptotic variance (4.5) 
with Vfc(a) from (4.6) instead of V k . Taking into account (4.16) or (4.17) we may replace 
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Vfc(a) by the reduced sum X^eG ^2/ ( a ) so ^ na ^ the secon d P ar t of (4.18) is a consequence 
of 



H E ^(«) ) 2 - °*(<0 <^ 2 §E (2n + 3)^ /^» _> . 

\ / Tr-"fe > fc— s-oo 



Here we have again used the notation Pj. and the standard covariance estimates from the 
very beginning of Step 4. Summarizing all Steps 1-4 completes the proof of Theorem 3.1. 

□ 



4.2 Proof of Lemma 3.1 

(2) (2) 

By definition of the signed measures 7^- and 7^ d in Section 2.2 and using algebraic induc- 
tion, for any bounded Borel-measurable function g : (R d ) 2 — > R 1 we obtain the relation 

A / [ g(x,y)^ ) red (dy)dx= [ g(x, y - x) 7 f \d(x, y)). (4.19) 



Let H + , H be a Hahn decomposition of M. d for Jx\ed> i- e -' 

7^(0 = ^Ze d (H + n (•)) and ^- d (.) = -4] red {H- n (•)) • 

We now apply (4.19) for g{x,y) = l Eo (x) 1 H +nE z {v) > where E » = i~h ^ + z for z G ^ ■ 

(2) 

Combining this with the definition of the (reduced) second factorial moment measures a x 
(2) 

and a x d of the unmarked PP X = X^i>i j see (2-5) for m = 2, and 

jP(A x B) = a'x ) (Ax B)- \ 2 \A\\B\ for all bounded A, B £ B(M. d ) , 

leads to 

X lx]red( H+ ^ E z) = J ^E {x)1 H+nEz {y-x)a^{d{x,y))-\ 2 \E \\H + ^E z \ 



= lE (X i )l H+nE .(X j - Xi) -EX(E )EX(H + C\E Z ). 

i,j>l 

Since o ^ H + n P 2 for z£Z d with |z| > 2 we may continue with 

X lx]red( H+ n = E^<5 Xj (£ ) X((H+ n £ z ) + X,) - EX(E )EX(H+ n £ z ) 

i>l 

= Ef(YX)-^f(YX) for |z|>2, (4.20) 

where 

f(Y, Y' z ) = S Xi (E Q ) X((H + n E z ) + Xi) < X(E ) X(E Z © E Q ) (4.21) 

i>l 

with F(-) = J2i>i frr, ((•) n E ) resp. *;'(•) = S Xj (((•) n £,) being restrictions 

of the stationary PP X = Y,i>i 6 Xi to Pq resp. E z ® E Q = [-1, l) d + z . Further, let Y and 
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Y z denote copies of the PPs Y and Y' z , respectively, which are assumed to be independent 
implying that Rf(Y,Y£) = EX(E )EX(H + n E z ) . Since Y is measurable w.r.t. A X (E ), 
whereas Y' z is ^4x(^ d \[ — (\ z \~ 1)> \ z \~ l] rf )-measurable, we are in a position to apply Lemma 2.2 
with/3(^4 x (£; o ),^4 x (R a! \[-(|z|-l),|z|-l] a! ) < /3 Xm (\z\ - §) for \z\ > 2. Hence, the estimate 
(2.10) together with (4.20) and (4.21) yields 

\^x,re d (H + nE z )\< 2^ x J\z\ - I)) 1 * 5 (maxjE/^^Fj) , Ef+*(Y,Y>) })*\ 

where the maximum term on the rhs has the finite upper bound 2 d ( 1+r i) E(X(E )) 2+2r > for 
5 = 2ij > in accordance with our assumptions. This is seen from (4.21) using the Cauchy- 
Schwarz inequality and the stationarity of X giving 

1 /2 

E/ 1+r >(Y,Y z ') < [E{X{E ) 2+2r ')E{X{[-l 1 l] d ) 2+2 ^ < 2 d( - 1+ ^ E(X (E ) 2+2r >) 

and the same upper bound for Ef 1+V (Y ,Y Z ) . By combining all the above estimates with 
X ^xled( H+ n [-§. D d ) ^ 3 d EX(E ) 2 we arrive at 

^ ( xU H+ ) * ^X(E ) 2 + 2^ (E(X(E o) r s )^ (^(M ~ \))^- 

z€Z d :\z\>2 

By the assumptions of Lemma 3.1 the moments and the series on the rhs are finite and the 

(2) 

same bound can be derived for — \j x red (H~) which shows the validity of (3.4). 

The proof of (3.5) resembles that of (3.4). First we extend the identity (4.19) to the (reduced) 
second factorial moment measure of the MPP Xm defined by (2.4) and (2.6) for m = 2 which 
reads as follows: 



A / J g(x,y)P° i x (C 1 xC 2 )a^ red (dy)dx = J g(x, y - x) P^(d x C 2 ) ag } (d(x, y)) 

R d R d (R d ) 2 

= E^fgiXi, Xj - X, ) l Cl (Mi) I Ca (Mj ) . 

i,j>l 

For the disjoint Borel sets G + and G~ defined by 

G + ^ = {xeR d : P^(dxC 2 ) > (^P^P^C,)} 

we replace g(x, y) in the above relation by g^(x, y) = 1e {x) TL E ±(y) , where E^ = G^C\E Z for 
\z\ > 2, and consider the restricted MPPs Y Q (-) = X M ((-) n (E Q x Ci)), Y£±(-) = X M {{-) n 
((Ef E Q ) x C 2 )) and their copies Y Q and Y' z ± , which are assumed to be stochastically 
independent. Further, in analogy to (4.21), define 

f(Y , Yi i± ) = ]T s (Xi,Mi)( E o x Ci) X M {(Et + Xi) x C 2 ) < X(E ) X(E Z E Q ) . 

It is rapidly seen that, for \z\ > 2 , 

Ef(Y X,±) = A J P° i x (C 1 xC 2 )a% ed (dx) and 

Ef(Y , Fj i± ) = EX M (E x d) EX A / x C 2 ) = A 2 P&(Ci) P&(C 2 ) |^ I 
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and in the same way as in the foregoing proof we find that, for \z\ > 2 

7,\ — 

2 



Ef(Y X,±)-^f(YoX,±)\< 2 d+1 {MX(E ) 2 + S )£s {f3 x J\z\ 



Finally, the decomposition ot x red {-) = ^ x red(') + -M ' I together with the previous estimate 
leads to 



A 

E, 



J |p^(d X C 2 ) - PUCl) P° M (C 2 )\ ^] red (dx) = mYoX,+) ~ ®f<XoX,+) 

Ef(Y Q , Y^_) - Ef(Y a X,-) ) " AP&(Ci) P° M {C2) U x ] red { E+ ) - 7 £ ed ( E z ) 



<2 d+2 (E(X(£ Q ) 2 + 5 ))^^ for \z\>2. 

Thus, the sum over all z 6 Z rf is finite in view of our assumptions and the above-proved 
relation (3.4) which completes the proof of Lemma 3.1. □ 

4.3 Proof of Theorem 3.3 

It suffices to show (3.6), since independent marks imply that P°j x (Ci x C 2 ) = P°j(Ci) P°j(C 2 ) 
for x 7^ o and any C±, C 2 G £>(M) so that the integrand on the rhs of (3.6) disappears which 
yields (3.7) for stationary independently MPPs. By the very definition of Y k {C) we obtain 
that 



Cav(Y k (C i ),Y k (C j )) = Ej2^w k (X P ){t Ci (M p ) - P^(C,)) (Ic 3 (M p ) - P^C 3 )) 

p>i 

+ W] E E **w k (X p )l Wk (X q ){l Ct (M p ) - P&{Ci)) (lo,(M g ) - P&CCj-)) • (4.22) 

p,q>l 

Expanding the difference terms in the parentheses leads to eight expressions which, up to 
constant factors, take either the form 

Ej2^W k (X p )l c (M p ) = X\W k \ PfciC) or E ^ *^W k {X p )l Wk {X q )l Ci (M p )l Cj {M q ) 

p>l p,q>l 



l Wh {x)lwMPT~ X (Ci x Cj)a^(d(x,y)) = \j P^iC, x Cj) lk {y) a% red (dy) , 

where y \— > "f k (y) = \ W k n (W k — y)\ denotes the set covariance function of W k . Summarizing 
all these terms gives 

Cov{Y k (C i )Xk(C j ))=x(P^(C i nC j )-P^(C i )P^(C j ))+X I ^ (ptfiCiXCj) 

R d 

- P^d) PtflCj x M) - P° M {C 3 ) P°?{G x M) + PUCi) P°m{C 3 )) a% ed {dx) . 
The integrand in the latter formula is dominated by the sum 

|P£f (Ci x Cj) - P^(Ci)P^(C 3 )\ + \PZf(Cj x M) - P^(C 3 )\ + \P°f(Ci x M) - , 

i 2 \ 

which, by (3.5), is integrable w.r.t. at red - Hence, (3.6) follows by (2.2) and Lebesgue's 
dominated convergence theorem. □ 
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4.4 Proof of Theorem 3.4 

We again expand the parentheses in the second term of the estimator (a : ^ 1 '*) fc defined by 
(3.8) and express the expectations in terms of P°f and offi j. Using the obvious relation 
7k{y) = / R d ^W k (x)^w k (y + x ) dx we find that, for any Ci,Cj £ B(M) , 

^ \(w k -x p )n(w k -x q )\ - J «x(d(^yJ) 

Pm i^i x Cj 



l Wk (x)l Wk (y + x) dxa® red (dy) = \f P^{d X Cj) a%j£y) • 



ik{y) 

R d R d R d 

As in the proof of Theorem 3.3 after summarizing all terms we obtain that 
n°S ] ) k = A(p^(a n Cj) - PUCi)P^iCj)) + A J (p°f{Ci x Cj) 

R d 

- P°f{C x M)P°j(Cj) - P^ x (Cj x M) J* (Q) + P&(C1) P^)) a^Cdx) , 

which, by comparing with (3.6), yields that E{a^) k = a^j . The asymptotic unbiasedness of 

(d% ] ) k is rapidly seen by the equality E{a£\ = Cov( Y k (Ci), Y k (Cj) ) = EY k (Ci)Y k (Cj) , 
which follows directly from (4.22), and (3.3). □ 

4.5 Proof of Theorem 3.5 

Since E(cr^ - (c?-^) fc ) 2 = Var(a^) fe + (oij - E(d^) k ) 2 we have to show that 

For notational ease, we put m(u,v) = (lc.(*u)-P^(Cj)) (l c ^ (u)-P°/(C,-)) , a fc = b k \W k \ 1/d , 

-) and r fc = r k (X p , X q ) m{M p , M q ) . 

lk{y~x) \ a k J ^ 

Hence, together with (2.3) and (3.1) we may rewrite (<7^ 3 ) fe as follows: 

&j\ = ^ n C i) + n Cj) - PUCi)P^{Cj)) + r k . (4.24) 

Using the definitions and relations (2.4) - (2.6) and J^ d r k (x,y + x)dx = w(\\y\\/a k ) we find 
that the expectation Er k can be expressed by 

(fxM) 2 R d M 2 
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The inner integral f M2 m(u, v) P°f (d(tt, v)j coincides with the integrand occurring in (3.6) 

(2) 

and this term is integrable w.r.t. a x > due to (3.5) which in turn is a consequence of (3.10) 
as shown in Lemma 3.1. Hence, by Condition (wb) and the dominated convergence theorem, 
we arrive at 

Et* — >A f f m(^)P^(d(u,«))ag red (d^ 



The definitions of A^ and Y k (-) by (2.3) and (3.1), respectively, reveal that EA/% = A and 
E Yfc(Cj Pi Cj) = . This combined with the last limit and (4.24) proves the first relation of 

(4.23) . To verify the second part of (4.23) we apply the Minkowski inequality to the rhs of 

(4.24) which yields the estimate 

(Var (d^ )J 1/2 < \Wu\- 1 ' 2 (Var Y k (Q n C,)) 1/2 + (Var X k ) 1/2 + (Var r fc ) 1/2 . 

The first summand on the rhs tends to as k — > oo since E Y k (C) 2 has a finite limit for any 
C € £>(M) as shown in Theorem 3.3 under condition (3.5). The second summand is easily 
seen to disappear as k — > oo if (3.4) is fulfilled, see e.g. [12], [16] or [17]. Condition (3.10) 
implies both (3.4) and (3.5), see Lemma 3.1. Therefore, it remains to show that Var — > 
as k — > oo . 

For this purpose we employ the variance formula (2.7) stated in Lemma 2.1 in the special 
case f(x,y,u,v) = rj~{x,y) m(u,v) . In this way we get the decomposition Varrfc = il + 
+ ij?' , where , and 1^ denote the three multiple integrals on the rhs of (2.7) 
with f(x, y, u, v) replaced by the product r k (x, y) m(u, v) . We will see that the integrals 

(2) (3) 

and IY are easy to estimate only by using (3.4) and (3.5) while in order to show that /} 
goes to as k — > oo, the full strength of the mixing condition (3.10) must be exhausted. 
Among others we use repeatedly the estimate 

for y€B(o,r w ), (4.25) 



lk{akV) \W k \ 

which follows directly from (2.2) and the choice of {b k } in (3.9). The definition of together 

(2) (2) 

with (4.25) and a X red (dx) = lx r ed^ x ) + Adx yields 

l4 X) | < 2 / {rk{^2)fay{d{ Xl ,x 2 ))=2\j^ )W 2 

m 2 w h ( Z ed m d ) + Xa d k f w 2 (||y||)dy) -^0, 



4A 

Wk. 



where the convergence results from Condition (wb) and (3.10), which implies iTxtedK^) < ~ 
oo by virtue of Lemma 3.1. Analogously, using besides (4.25) and Condition (wb) the relations 

<mwK[-\a k r w ],\a k r w ]]4 x ) and W k Q [jE z with H k = H k UdH k , 

z£H k 
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with the notation introduced at the beginning of the proof of Theorem 3.1, we obtain that 
|^ 2) | ^ 4 / r k (xi,x 2 ) r k (xi,x 3 ) {d(x 1 ,x 2 ,x 3 )) 



,/ys 



■') 

16: 

< 



16 771^ 

j^nf E a l((^ © [-\a k r w ],\a k r w ]] d ) x (£ 2 [-\a k r w ],\a k r w ]] d ) x £ 2 ) . 



Since the cube E z © [— [a^r^], [afcr w ]] d decomposes into (2\a k r w ~\ + l) d disjoint unit cubes 
i%\E zl xE Z2 xE Z3 ) 



and a^y {E Z1 x E Z2 x ^ 3 ) < K(X(E )) 3 by Holder's inequality, we may proceed with 



Here we have used the moment condition in (3.10), (4.3), and the assumptions (3.9) imposed 
on the sequence {b k } . 

In order to prove that /[ vanishes as k — > oo we first evaluate the inner integrals over the 
product m(u±, u 2 ) m(u 3 , ti 4 ) with m(u,v) = (ic^it) — -Pj^(Cj)) (la, ) — P°i(Cj)) so that 
1^ can be written as linear combination of 16 integrals taking the form 

Jk = / r k (xi,x 2 )r k (x 3 ,x 4: ) P^' X2 ' X3 ' X4 ( x D^a^ (d(x 1 ,x 2 ,x 3 ,x 4 )) 



r=l 



P^ 2 (Z?i x D 2 )P X ?> X \D 3 x J D 4 )ag ) (d(x 1 ,x 2 )) (d(x 3 , x 4 )) 



4 r fc (xi,x 2 )r fc (x3,x 4 ) {orx M ~ a x M x axL)( d ( a;i '' ltl '-' a;4 '' lt4 )) ' 

x (R d xD r ) 

r — l 



where the mark sets D\,D 3 G {Cj,M} and D 2 ,D 4 G {Cj,M} are fixed in what follows and 

(4) _ (2) (2) 

■X-M -X-M X-M 



the signed measure at^ — x on i3((M d x M) 4 ) (and its total variation measure 



at*x ~ a x' ' x I ) come m t° pl & y by virtue of the definition (2.6) for the ?n-point Palm 



mark distribution in case m = 2 and m = 4 . 



As |zi — z 2 \ > \a k r w ~\ (where, as above, \z\ denotes the maximum norm of z G Z d ) implies 
11^2 — ^i|| > o, k r w and thus r^(xi,X2) = for all x\ G E Zl ,x 2 G E Z2 , we deduce from (4.25) 
together with Condition («;&) that 



2 / 2 K^1 \ 

E + E E E y — 4 , (4.26) 

\ n=0 n>2Krv,l / (ai,« 2 )eSfc (.23 ,24) £5* 



where S*,. = {(u,v) G H k xH k : \u-v\ < \a k r w ]} , Sfc^z) = {(zi,z 2 ) G S fc : mm = n} 

and y zii22i23)24 = ~ a x M x "xL I ( x r=l( £7 2r x A-)) for any Zi,...,z 4 G Z d . 

Obviously, for any fixed 2 G iifc, at most 2(2 [a^r,,,] + l) d (4 [a^ru,] + l) d pairs (z 3 , z 4 ) be- 
long to (Jnl^o Sk,n{z) and the number of pairs (21,^2) m S k does not exceed the product 
j£H k (2 [a^r^] + l) d . Finally, remembering that a k = b k {W^ 1 ^ and using the evident 
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estimate V Z1 z% Z3 Z4 < 2 E,(X(E )) 4: together with (4.3) and Condition (wb), we arrive at 

U E ! E E v mu <c m *^(4m)\^j. 

It remains to estimate the sums on the rhs of (4.26) running over n > 2 [a^r^] . For the signed 
measure a^| f — x we consider the Hahn decomposition H + ,H~ G B((M. d x M) 4 ) 

yielding positive (negative) values on subsets of H + (H~). Recall that K a = [— a, a] d . For 
fixed z\ G H^, 22 G Hk H (Kr akr i + z\) and (2:3,2:4) G Sk >n (^l), we now consider the 
decompsition V^,^ = ^ + V^,^ with 



= ± («£, " «g, x 41) * n _i (A, x AO) • 



r=l 



Since (23,24) G Sk, n (zi) means that 23 G iffc Pi (K^ + 21), where K c a = R \ i"T a , and 
2:4 G fl (i^|- aj . rm -| + 23) n (if^ + 21) , we define MPPs Yfc and Y^ as the restrictions of Xm 
to (K^ ak r w ~\+l/2 + z i) x M and (-^_]/ 2 + z i) x respectively. Let furthermore Yfc and Y^ 
be copies of Y& and Y^ which are independent. 

Next we define functions / + (Yfc, Y^) and / _ (Yjfc, Y^) by 

/ ± (Y fc , Y r = Yf ^ft' M P' X ^ M ^ X 's,M' s ,X' tl M[) , 

p,q>l s,t>l 

, 4 

where I±(- • • ) denote the indicator functions of the sets H fl x (A r x A) so that we get 

r=l 

^f, Z2 ,*3,*4 = E/ ± (Y fc , 1*) - E/ ± (Y fc , Y£ for (2 l5 2 2 ) G 5 fc , (23, 24) G S k , n ( Zl ) . 



Hence, having in mind the stationarity of Xm, we are in a position to apply the covariance 
inequality (2.10), which provides for 77 > and n > 2 [afer w ] that 

V 

r± nl at _ I ~ N AtTSC 1 _ \\ \ !+»> 



< 2f /3(^(K K , wl+ i/ 2 + Zx),A(K^_ 1/2 + 21)) 
2 4 
E( ]J *m(A. x D r )) 2+2n E( XM{E Zr x AO) 



2+2r? "Q v 
r=l r=3 

< 2(/3^>- K^] - 1))^ (E(X(A,)) 4 + 4 ")^ . ( 4 .27) 



In the last step we have used the Cauchy-Schwarz inequality and the definition (2.9) of the 
/3-mixing rate Px M ' Finally, setting r] = 5/4 with 5 > from (3.10) the estimate (4.27) 
enables us to derive the following bound of that part on the rhs of (4.26) connected with the 
series over n > 2 [a^r^] : 

c n **** (2\a k r w ] + lf d H 2n + l ) d ~ ( 2n ~ !)") (P*x M (n ~ \a k r w ] - 1))^ . 

n>2\a k r nl ~] 

Combining a k = 6fc|Wfc| 1 / <i an d (4.3) with condition (3.10) and the choice of {b k } in (3.9), it 
is easily checked that the latter expression and thus Jj, tend to as k — > 00 . This completes 
the proof of Theorem 3.5. □ 
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5 Examples 

5.1 m-dependent marked point processes 

A stationary MPP Xm is called m-dependent if, for any B,B' £ B(M. d ), the a— algebras 
Ax M (B) and Ax M (B') are stochastically independent if inf{|x — y\ : x £ B, y £ B'} > m 
or, equivalently, 

f3(Ax M {K a ),Ax M {Ka+b))=Q for fo > m and a>0. 

In terms of the corresponding mixing rates this means that (3 X (r) = /3 X (r) = if r > 
m. For m-dependent MPPs Xm it is evident that Condition (5(8) in Theorem 3.1 is only 
meaningful for 5 = 0, that is, E(X([0, l] rf )) 2 < oo . This condition also implies (3.4) and 
(3.5). Likewise, the assumption (3.10) of Theorem 3.5 reduces to E(X([0, l] d )) 4 < oo which 
suffices to prove the Inconsistency of the empirical covariance matrix sj^ . 



5.2 Geostatistically marked point processes 

Let X = ^2 n >iSx n be an unmarked simple PP on M. d and M = {M(x), x £ M. d } be a 
measurable random field on M. taking values in the Polish mark space M. Further assume that 
X and M are stochastically independent over a common probability space (O, A, P). An MPP 
Xm — Y2 n >l fi(X n ,M n ) with atoms X n of X and marks M n = M(X n ) is called geostatistically 
marked. Equivalently, the random counting measure Xm £ Nm can be represented by means 
of the Borel sets Af -1 (C) = {x £ R d : M{x) £ C} (if C £ B(M)) by 

X M (BxC) = X(BnM- 1 {C)) for £ x C £ B(R d ) x B(M) . (5.1) 

Obviously, if both the PP X and the mark field M are stationary then so is Xm and vice 
versa. Furthermore, the m-dimensional distributions of M coincide with the m-point Palm 
mark distributions of Xm- The following Lemma allows to estimate the /3-mixing coefficient 
(2.8) by the sum of the corresponding coefficients of the PP X and the mark field M. 

Lemma 5.1. Let the MPP Xm be defined by (5.1) with an unmarked PP and a random mark 
field M being stochastically independent of each other. Then, for any B, B' £ B(M. d ) , 

P(A Xm (B),A Xm {B')) < p(Ax(B),A x (B')) + /3(A M (B),A M (B')) , (5.2) 

where the a— algebras Ax(B), Ax(B') and Am (B), Am (B 1 ) are generated by the restriction 
of X and M , respectively, to the sets B ,B'. 

To sketch a proof for (5.2), we regard the differences A(A i ,A' j ) = P(A; n Aj) - P(Ai) P(A' j ) 
for two finite partitions {Ai} and {Aj} of Q consisting of events of the form 

k e 
Ai = f]{X M (B p x C p ) £ r P)i } , A'j = f]{X M (B' q x C' q ) £ T' qj } with r^,^ c %\ , 

p=l q=l 

with pairwise disjoint bounded Borel sets B\, B^ Q B and B[, B^ C B' . Making use of 
(5.1) combined with the independence assumption yields the identity 

A(^, 4) = J J (v Ax{mAx[BI) - F Ax{B) x F Ax{BI) yAi n A'j) dP AM{mAM{BI) 
n n 

+ J J F Ax(B)i A i) F Ax(B>){A'j) d(^Au{B)9Au(B') ~ F A M (B) X F Am(B')) , 

a n 
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which by (2.8) and the integral form of the total variation confirms (5.2). 
5.3 Cox processes on the boundary of Boolean models 

Let H = Un>i(^™ + ^n) De a Boolean model, see e.g. [20], governed by some stationary 
Poisson process Xm>i $Y„ m ^ with intensity A > and a sequence {H n } n >i of independent 
copies of some random convex, compact set So called typical grain (where we may assume that 
o € Ho). With the radius functional ||So|| = sup{||a;|| : x 6 Ho}, the condition E||Ho|| d < oo 
ensures that H is a well-defined random closed set. We consider a marked Cox process Xm, 
where the unmarked Cox process X = Y2n>i ^ x n is concentrated on the boundary <9H of 
H with random intensity measure being proportional to the (d — l)-dimensional Hausdorff 
measure 7id-i on 9E,. As marks M n we take the outer unit normal vectors at the points 
X n 6 83, which are (a.s.) well-defined for n > 1 due to the assumed convexity of Ho- This 
example with marks given by the orientation of outer normals in random boundary points 
may occur rather specific. However, this way our asymptotic results may be used to construct 
asymptotic tests for the fit of a Boolean model to a given dataset w.r.t. its rose of directions. 
For instance, if the typical grain Ho is rotation-invariant (implying the isotropy of H), then 
the Palm mark distribution P M of the stationary MPP Xm = Y2 n >l ^(X n ,M n ) is the uniform 
distribution on the unit sphere S d_1 in R d . We will now discuss assumptions ensuring that 
Condition (3(5) and (3.10) hold, which are required for our CLT (3.2) and the consistent 
estimation of the covariances (3.3), respectively. Using slight modifications of the proofs for 
Lemmas 5.1 and 5.2 in [15] one can show that for a, b > 

(3(Ax M (K a ),Ax M (K c a+b ))<\2 d+2 ( 3 + y) d ^(llHoflOHoll^}). 

According to (2.9) we may thus define the /3-mixing rates /3 XM ( r ) an d Px M ( r ) ^ or r — I ^y 
0* x Jr) = c 12 E(||Ho|| d I{||S ||>r/4})> sup f3(Ax M (K a ),Ax M (K c a+r )), 

l/2<a<r 

Fxjr) = -^EOlSortllHoll > r/4}) > supa"^ 1 ) P{Ax M (K a ),Ax M (K c a+r )) , 

' a>r 

where c 12 = X2 d+2 l d ~ l . 

It is easily seen that E||H || M < oo implies r 2d ~ x ^ (r) — > 0. Moreover, E||H || 2d(p+ ^/ <5+E < 

oo for some e > ensures r d ~ l (/3 Xm (r)) 5 ^ 2p+<5 ^dr < oo for p > 0. Since the random 
intensity measure of X on E a and thus also X(E a ) has moments of any order by virtue of 
E||Ho|| d < oo , the parameter 5 > in Condition (3(5) and (3.10) can be chosen arbitrarily 
large. This results in the following lemma. 

Lemma 5.2. For the above-defined stationary marked Cox process Xm on the boundary of 
a Boolean model H with typical grain Ho the assumptions of the Theorems 3.1 and 3.5 are 
satisfied whenever 

E||H || M+£ < oo for some e > . (5.3) 

Remark: The marked Cox process Xm is m-dependent if ||Ho|| is bounded by some constant. 
By using approximation techniques with truncated grains as suggested in [15], pp. 299-302, 
it can be shown that (5.3) is just needed for e = . Moreover, the statistical analysis of roses 
of directions via marked Cox processes applies also in case of non-Boolean /3-mixing fibre 
processes, see e.g. [12] for Voronoi tessellations. 
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6 Simulation study 

6.1 Moving average model in IR 2 

In this section we introduce an m-dependent MPP model, which was used for our simulations 
since it allows to control the range of spatial dependence for a fixed Palm mark distribution. 
The locations of this MPP are given by a homogeneous Poisson process Xm>i ^X n m Each 
point is marked by a direction in the upper half of the unit circle. In order to construct 
the Palm mark distribution, we consider the projected normal distribution PN2(a, n) on 
S 1 . By definition, Y ~ PN2(a, k) means that Y = -v^n for some Gaussian random vector 

Z ~ N~2{ a J K ) m H^ 2 with an invertible covariance matrix k. Note that PN2(o, cr 2 ^) is 
the uniform distribution on S 1 for all a 2 > 0, where I2 is the identity matrix. Formulas 
for the density of a projected normal distribution can be found in [19]. Let {M^ } n >i 
be iid N2(o, /^-distributed random vectors. The stability of the normal distribution w.r.t. 
convolution yields 



Mi 2 ) 



ZZiM?\ llXi - Xnl \< p} _.p N2(0jK)) 



E£i M i 1):l[ {ii^-x n ||<p} 



for any p > controlling the range of dependence. The marks of our model are finally defined 
as the axial versions M n = M n 1§+ {M n 2 ^ ) — M n !§- (M n ) of the averages M n , i.e. , points 

on the lower half-circle §1 are rotated by tt. Due to the moving average approach defining the 

(2) 

preliminary marks {M n }, we call the MPP Xm = Sn>i ^(X„,M„) the moving average model 
(MAM). The MAM is clearly m-dependent, where the range of dependence is controlled by 
the averaging parameter p. 



6.2 Tests 

By simulations of the MAM we investigated the performance of the asymptotic x 2 -g°°dness- 
of-fit test, which is based on the test statistic 

Tk = Y I^k lY k Ke- 
lt— >oo 

If is chosen as the L 2 -consistent estimator (a^ and (-Pj^)o denotes a hypothetical 

Palm mark distribution, the hypothesis Hq : P°j = (-Pj^)o is rejected, if > xj i_ a > where a 
is the level of significance, and x| , denotes the 1 — a-quantile of the ^-distribution. This 
test will be referred to as 'test for the typical mark distribution' (TMD). The construction of 

(a^j ) k involves the sequence of bandwidths We set 

3 

bk = c|Wfc|~4d for a constant c > 0, (6.1) 

and chose c such that condition (wb) for the inconsistency of (3y ) fc was satisfied. The 
specific choice of the constant c is crucial for test performance, as discussed below. This choice 
of c can be avoided if S is not estimated from the data to be tested but incorporated into Hq. 
This means, we specify an MPP as null model, such that So is either theoretically known 
or otherwise can be approximated by Monte-Carlo simulation. By means of the combined 
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null hypothesis Hq : P M = {Pm)q an d ^ = ^o> the test exploits not only information on 
the distribution of the typical mark but additionally considers asymptotic effects of spatial 
dependence. The test can thus be used to investigate if a given point pattern differs from 
the MPP null model w.r.t. the Palm mark distribution. We will therefore refer to it as 'test 
for mark-oriented goodness of model fit' (MGM). By the strong law of large numbers and the 

^ to) 

asymptotic unbiasedness of (a^ ) k , a strongly consistent Monte-Carlo estimator for "Sq in 
an MPP model Xm is given by 



where X M , . . . ,Xj^ are independent realizations of Xm- Thus, for large k and n the test 
statistic Tfc jn = Y^S^^Y^ has an approximate xf distribution. If a is the level of signif- 
icance, the MGM test rejects Hq, if T^ n > Xei-a- ^ ne estimator ^k,n can also be used 
to construct a test for the typical mark distribution if independent replications of a point 
patterns are to be tested. In that case X M \ ...jXj^ are the replications. Note that for 
replicated point patterns, Hq does not incorporate an assumption on S and hence the cor- 
responding test differs from the MGM test. The edge-corrected unbiased estimator 

was not used for the Monte-Carlo estimates in our simulation study, since (a^ ) fc can be 
computed more efficiently 



6.3 Model parameters 

The MAM was simulated on the observation window — [—1500, 1500] 2 . The expected 

number of points was set to EA(Wi5oo) = 3125. The asymptotic behavior of the test was 
studied by considering smaller observation windows corresponding to an expected number 
of 300, 600, . . . , 3000 points. Spatial stochastic dependence of marks was varied by the 
parameter p G {0,50, . . . ,300}. In the MAM, marks of points with distance no larger than 
2p in general exhibit stochastic dependence. If, on the contrary, two points are separated 
by more than 2p, their marks are independent. Thus, p = corresponds to independent 
marking. Deviations of the projected normal distribution from the uniform distribution 
on were controlled by varying K\2 £ {0,0.1,0.2,0.4,0.8}, where K\2 = represents the 
uniformly distributed case. The parameter k\\ = K22 = 1 was kept constant. The bins 
Ci, . . . ,Ce G £>(§+) for the x 2 -goodness-of-fit test were chosen as 



We will discuss the case I = 8, where the bins had a width of 20°. Simulations for i = 17 did 
not reveal different general effects. 

6.4 Simulation results 

All simulation results are based on 10000 model realizations per scenario. Type II errors were 
computed for realizations where K12 7^ 0, which means that the mark distribution was not 
uniform on whereas Hq : P°j = U(B>\) hypothesized a uniform Palm mark distribution 
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on §^ (corresponding to K\2 = 0). 

The performance of the MGM test is visualized in Tab. 1. Empirical type I errors of the 
MGM test were close to the theoretical levels of significance for a = 0.025,0.05, and 0.1 
with maximum deviations of around 0.015. They were hardly affected by the dependence 
parameter p. Type II errors increased with p. Under independent marking (p = 0) as well 
as for p = 50, error levels were close to for K\2 € {0.2,0.4,0.8}. However, for an extreme 
range of dependence (p = 300) even for a strong deviation of the data from a uniform Palm 
mark distribution (k±2 = 0.4), rejection rates were only between 30 and 40%. For p = 300 
the range of dependence corresponds to 20% of the sidelength of W. 

Experiments with the TMD test revealed that the choice of the bandwidth parameter c in 
(6.1) is critical for test performance (Tab. 3). Whereas large values of c result in small type 
I errors, they decrease the power of the test. On the other hand, small values for c lead 
to superior power but increase type I errors (Tab. 3). The empirical errors in Fig. 2 were 
computed for c = 50 which yielded a reasonable compromise with respect to the two error 
types. In comparison to the MGM test the TMD test exhibits a higher sensitivity of empirical 
type I errors for varying values of K\2, i.e., w.r.t. deviations from the uniform distribution 
on §1. Moreover, type II errors of the TMD test were up to 20% higher than for the MGM 
test. 

Tab. 3 and Fig. 1 illustrate test performance w.r.t. the mean number of points in W and 
the dependence parameter p. The simulation experiments were conducted for a = 0.05. 
For power analysis, the tested data was simulated for K12 = 0.4, and thus the Palm mark 
distribution strongly differed from the uniform distribution on S 1 of Hq. At a mean number of 
3000 observed points, Hq was reliably rejected by the TMD test once p < 150 (for c = 50). For 
p < 100 already 2000 expected points were sufficient to reject Hq for almost all realizations. 
The MGM test required around 500 points less than the TMD test in order to achieve 
comparable rejection rates (Fig. 1). 

In summary, our simulation results indicate that the MGM test outperforms the TMD test 
especially with respect to power. This result is plausible since the additional information 
incorporated into Hq by specification of a model covariance matrix can be expected to result 
in a more specific test. It seems difficult to derive a general rule of thumb relating the required 
size of the observation window to the dependence structure of the data and the intensity of 
the point pattern. However, Fig. 1 and Tab. 3 provide an idea on the practical requirements 
for asymptotic testing. 
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Figure 1: Empirical errors of types I and II for the MGM test plotted against the mean 
number of points in the observation window (K12 = 0.4, a = 0.5). Different colors correspond 
to different values of the dependence parameter p. 
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Type II error 
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.0.1 
0.2 
0.4 
0.8 



0.025 0.05 



0.1 



Kf2 

o 
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0.8 



0.025 0.05 



0.1 



50 



Q. "* 

>■ O 



K12 
o0.1 
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1^2 

> 
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p = 300 



K12 
o0.1 

0.2 
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0.025 0.05 



0.1 



Table 1: Empirical errors of types I and II for the MGM test. 
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Type I error 



Type II error 



Kj 2 

■ 
0.1 
0.2 
0.4 
0.8 



0.025 0.05 



0.125 



0J 



K, 2 
o0.1 

0.2 
+ 0.4 
0.8 



0.025 0.05 



0.125 



K^2 

. 
0.1 
0.2 
0.4 
0.8 



0.025 0.05 



0.125 



50 



Kj 2 

= 0.1 
0.2 

h0.4 
0.8 



0.025 0.05 



0.125 



K<2 

■ 
0.1 
0.2 

■ 0.4 
0.8 



0.025 0.05 



0.125 
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Table 2: Empirical errors of types I and II for the TMD test (c = 50). 
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Table 3: Empirical errors of types I and II for the TMD test plotted against the mean 
number of points in the observation window (K42 = 0.4, £ = 8, and a = 0.05). Different 
colors correspond to different values of the dependence parameter p. 



33 



